NASA/CR-2000-209853 
ICASE  Report  No.  2000-5 


Formulation  of  a  Two-scale  Model  of  Turbulence 


Robert  Rubinstein 
ICASE,  Hampton,  Virginia 


Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  VA 

Operated  by  Universities  Space  Research  Association 


National  Aeronautics  and 
Space  Administration 


20000308  049 


Langley  Research  Center  Prepared  for  Langiey  Research  Center 

Hampton,  Virginia  23681-2199  under  Contract  NAS  1-97046 

February  2000 


distribution  statement  a 

Approved  for  Public  Release 
Distribution  Unlimited 

DTIC  QUALITY  INSPECTED  5 


FORMULATION  OF  A  TWO-SCALE  MODEL  OF  TURBULENCE 

ROBERT  RUBINSTEIN* 

Abstract.  A  two-scale  turbulence  model  is  derived  by  averaging  the  two-point  spectral  evolution  equa¬ 
tion.  In  this  model,  the  inertial  range  energy  transfer  and  the  dissipation  rate  can  be  unequal.  The  model 
is  shown  to  reduce  to  a  standard  two-equation  model  in  decaying  turbulence. 
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1.  Introduction.  Despite  the  remarkable  success  of  the  two-equation  turbulence  model  in  predicting 
many  practically  important  turbulent  flows,  it  has  some  important  shortcomings.  For  example,  the  linear 
eddy-viscosity  relation  between  the  Reynolds  stress  and  the  mean  strain  rate  is  inadequate  in  many  problems; 
nonlinear  eddy  viscosity  models,  algebraic  models,  and  finally  the  Reynolds  stress  transport  model  [1],  have 
been  developed  in  response  to  this  problem. 

Whereas  these  models  focus  their  attention  on  the  Reynolds  stresses  and  leave  the  two-equation  model 
itself  basically  intact,  a  complementary  line  of  research  [2],  [3]  has  attempted  to  improve  the  two-equation 
model  by  addressing  the  over-simplification  inherent  in  any  description  of  turbulence  by  a  single  length- 
scale.  It  lead  to  the  formulation  of  multiple  scale  models  of  turbulence  in  which  the  transport  equations 
for  turbulence  kinetic  energy  and  dissipation  rate  are  each  replaced  by  transport  equations  for  the  kinetic 
energy  and  dissipation  rate  pertaining  to  a  definite  range  of  scales  of  motion. 

Multiple  scale  modeling  attempts  to  treat  the  response  of  turbulence  to  changes  in  the  large-scale  motion 
more  realistically  than  the  two-equation  model.  Whereas  the  two-equation  model  assumes  that  the  inertial 
range  can  adjust  instantaneously  to  changes  at  the  large  scales,  multiple  scale  modeling  allows  time-delays 
in  this  response  and  thereby  permits  a  more  refined  picture  of  the  time  dependence  of  turbulence. 

A  representative  multiple  scale  model  is  the  model  of  Hanjalic,  Launder,  and  Schiestel  (HLS,  [2]),  which 
takes  the  form  for  homogeneous  turbulence 

kp  —  P  ep 
kt  —  tp  —  et 

ep  =  CplefP-CpS 

•^p  ftp 

(1-1)  it  =  Ctl¥±-Ct  S 

kt  kt 

In  this  model,  the  turbulent  fluctuations  are  partitioned  into  two  regions  identified  by  the  subscripts  p 
(production)  and  t  (turbulence)  of  large-  and  small-scale  fluctuations  respectively.  Otherwise,  the  standard 
notation  is  used  in  Eq.  (1.1):  k  denotes  turbulence  kinetic  energy,  e  is  the  dissipation  rate,  P  is  production, 
and  Cpi ,  Cp2 ,  Cti ,  Ct2  are  model  constants.  Modeling  the  behavior  of  different  scales  of  motion  abandons 
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the  strictly  single-point  approach;  it  attempts  a  compromise  between  two-point  models,  which  offer  greater 
physical  fidelity  at  the  expense  of  greater  computational  complexity,  and  single-point  modeling. 

Early  research  on  this  type  of  model  was  frustrated  by  the  inevitable  appearance  of  a  large  number  of 
model  constants  for  which  calibrating  experiments  could  not  be  readily  identified;  thus,  Eq.  (1.1)  requires 
four  constants  in  place  of  the  two  constants  of  the  comparable  standard  two-equation  model.  However, 
interest  in  multiple  scale  models  has  recently  been  revived  by  the  possibility  that  numerical  simulations 
could  help  identify  the  model  constants  [4].  At  the  same  time,  new  approaches  to  multiple  scale  modeling 
based  on  two-point  models  have  been  advanced  [3],  [5],  [6]  which  promise  to  eliminate,  or  at  least  reduce, 
the  required  empirical  input. 

The  goal  of  the  present  work  is  to  derive  a  two-scale  model  in  which  energy  transfer  and  dissipation 
can  be  distinct  and  satisfy  different  transport  equations.  Kolmogorov’s  principle  of  locality  of  inertial  range 
transfer  [7]  is  used  to  separate  these  effects.  Unlike  the  previous  model  of  this  type  [6],  the  coefficients  in 
the  present  model  are  independent  of  Reynolds  number. 

The  multiple-scale  viewpoint  permits  a  fresh  derivation  of  the  two-equation  model  [3],  [5],  [6].  The 
two-equation  model  will  be  reconsidered  from  this  viewpoint.  Then  a  two-scale  model  is  derived  in  which 
transfer  through  the  large  scales  and  the  dissipation  rate  can  be  unequal. 

Multiple-scale  effects  are  expected  to  be  important  in  turbulent  flows  dominated  by  disequilibrium 
between  large  and  small  scales.  When  such  effects  are  absent,  as  they  are  in  the  standard  self-similar  flows 
used  to  calibrate  turbulence  models,  the  multiple-scale  must  reduce  to  a  single-scale  model  with  a  unique 
similarity  solution.  This  principle  is  applied  to  both  the  HLS  model  and  the  proposed  model  in  the  case 
of  decaying  turbulence.  Conditions  are  found  which  prevent  the  existence  of  more  than  one  exponent  for 
power-law  decay.  In  both  cases,  the  condition  is  a  simple  inequality  among  the  model  constants. 

The  occurence  of  partitioned  values  of  turbulence  kinetic  energy  and  dissipation  rate  naturally  generates 
many  new  possibilities  in  relating  the  Reynolds  stress  and  the  strain  rate;  however,  this  issue  will  not  be 
addressed  in  the  present  work. 

2.  The  Two-equation  Model  as  a  Single-scale  Model.  The  starting  point  is  the  isotropic  part  of 
the  spectral  evolution  equation  [9] 

(2.1)  E{k)  =  T(k)  +  n(«)  -  D(k)  +  V{k) 

where  E(k)  is  the  energy  spectrum,  T(k)  is  energy  transfer  spectrum,  U(k)  is  the  production  spectrum, 
D(k)  is  the  dissipation  spectrum,  and  V(k)  represents  diffusion  effects.  Only  homogeneous  turbulence  will 
be  considered,  so  that  V(k)  =  0.  The  dissipation  spectrum  is  defined  by 

(2.2)  D(k)  =  vk2E(k) 


and  closures  must  be  provided  for  T(k)  and  U(k). 
The  integrated  quantities  are 


(2.3) 


E{k)(Ik  = 
T(k)cIk  — 
n  (tt)dft  = 
D(k)cIk  = 


k 

0 

P 
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The  second  relation  expresses  the  conservation  of  energy  by  nonlinear  interaction.  The  energy  equation 
follows  from  integrating  Eq.  (2.1)  over  all  wavenumbers.  Using  Eqs.  (2.3), 

(2-4)  k  =  P-e 

This  result  is  independent  of  the  analytical  form  of  E(k)  and  II (/c),  and  of  the  closure  used  to  define  T(k). 

To  find  a  second  equation,  it  will  be  necessary  to  introduce  some  specific  assumptions.  The  simplest 
steady  solution  of  Eq.  (2.1)  is  the  Kolmogorov  steady  state  defined  by 


(2.5) 


corresponding  to 


E(k)  =  ^ 


0  if  K  <  Ko 

CkT2^ ft-5/3  if  Kq  <  K  <  Kd 


l  o 


if  K  >  Kd 


(2-6)  T(k)  =  -T{6(k  -  «o)  -  S(k  -  Kd)} 

where  T  represents  the  energy  transfer  through  the  inertial  range.  In  Eqs.  (2.5)-(2.6),  k0  is  the  integral 
scale,  and  Kd  ~  (e/V3)1/4  is  the  Kolmogorov  scale.  To  maintain  a  steady  state,  production  must  balance 
transfer  into  the  inertial  range,  and  transfer  out  of  the  inertial  range  must  balance  viscous  dissipation,  so 
that 


(2.7)  P  =  e  =  T 

Moie  complex  analytical  forms  for  E(k)  could  be  introduced,  but  this  change  adds  nothing  essential. 

Now  generalize  Eq.  (2.5)  to  the  time-dependent  case  by  letting  e  =  e{t)  and  k0  =  K0(t ).  For  high 
Reynolds  number  turbulence,  ignore  the  consequent  evolution  of  n ^  can  be  ignored.  Assume  that  the  second 
equality  from  Eq.  (2.7),  namely 

(2-8)  e  =  T 

continues  to  apply.  Since 

(2-9)  k=\cKe^K-^ 

it  follows  that 


(2.10) 


k  =  CK{t-1/3Kf/3i  -  e2/3«0-6/3«o) 


but  this  equation  does  not  lead  to  the  desired  second  equation  directly,  because  it  contatins  the  new  unknown 
k0.  We  consider  two  methods  by  which  a  second  equation  can  be  derived. 


2.1.  The  method  of  Schiestel.  The  problem  is  solved  in  [3]  by  postulating  that 


(2.11) 


Kq  OC 


E(k  o) 


which  is  equivalent,  for  a  Kolmogorov  spectrum,  to 
(2-12)  kq/aco  =  'ye/ k 
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Then  Eqs.  (2.10)-(2.12)  give 

(2.13) 

which  can  be  re-arranged  as 

(2.14)  «  = 

with  a  rather  good  value  for  C€ i  and  a  value  for  C€ 2  which  depends  on  the  choice  of  7.  This  is  essentially 
Eq.  (27)  of  [3]. 

2.2.  The  method  of  moments.  It  can  reasonably  be  objected  that  by  assuming  the  length-scale 
equation  Eq.  (2.12)  the  problem  has  simply  been  transferred  to  another  variable  and  therefore  Eq.  (2.14) 
has  not  been  derived,  but  simply  postulated  indirectly. 

An  interesting  alternative  has  been  suggested  in  [5],  [6]:  by  taking  another  moment  of  the  spectral 
evolution  equation  Eq.  (2.1),  an  additional  equation  is  obtained  and  ko  can  be  eliminated  between  this 
equation  and  Eq.  (2.10).  This  implies  a  reduction  of  Eq.  (2.1)  to  a  finite-dimensional  system  by  a  Galerkin 
approximation. 

For  example,  multiply  Eq.  (2.1)  by  k~ 1  and  integrate.  With  the  assumptions  Eqs.  (2.5)— (2.6) , 


±\f 

dt  J0 


00 


=  o5/3 


(2.15) 
and 

(2.16) 
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—  ^  Ck  € 

Ko 

5 

k 

2  k  . 

K0 

5  ck0 

1/3  Jf-5/3, 


f 


}  _  T(k)  6  €  € 

dK  — —  = +  —  » 


K  K  0  Kd  « 0 

To  make  the  calculation  definite,  assume  for  the  production  spectrum 


(2.17) 


(  0 


|  CDe'/3K-V3S2 
0 


if  K  <  Kq 
if  Kq  <  K  <  Kd 
if  K  >  Kd 


where  S  is  the  mean  strain  rate.  Eq.  (2.17)  is  consistent  with  the  usual  two-equation  model,  since  integration 
over  k  leads  to 

(2.18)  P  =  Coe^K-^S2  =  ^%—S2  =  C„—S2 

9  C  e  e 

which,  for  the  values  of  inertial  range  constants  recommended  by  Yakhot  and  Orszag  [8]  Cd  ~  0.5,  Ck  ~  1-6, 
gives 


(2.19) 


i  _  4  Cd  ^  4 
”  9  Cj£  ~  45 


With  this  choice, 

(2.20) 


n(«) 

K 


=  -7Cn^W/3^ 


4  P_ 

7  Ko 
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Note  that 


(221) 


f 


D(k)  2/3 

(IK  - -  ~  VK^ 


'  Re~1^2 


is  negligible  in  the  high  Reynolds  number  limit. 

Combining  the  results  Eqs.  (2.15),  (2.16),  (2.20),  and  (2.21), 


(2.22) 


«o  5  e 


_L  i  — 

K0  7  Ko 


consequently, 

(2-23)  -  \W 

with  a  value  Ce\  =  15/14  which  is  somewhat  too  small,  but  also  Ce 2  =  0. 

The  absence  of  a  destruction  term  comes  about  because  of  the  cancellation  of  k  and  e  in  Eq.  (2.22). 
This  cancellation  proves  to  be  independent  of  the  order  of  the  moment  taken,  provided  of  course  that  the 
integrals  converge  at  large  k. 

The  question  then  arises  whether  the  result  of  Eq.  (2.16)  might  not  depend  strongly  on  the  assumptions 
made  and  whether  more  generally,  we  should  not  obtain 

(2.24)  r  dK^-  =  -CT  — 

Jo  K  &0 

with  Ct  7^  1.  But  Eq.  (2.16)  is  reproduced  exactly  if  the  closure  for  T(k)  given  by  Eq.  (2.6)  is  replaced  by 
Leith  diffusion  closure  for  T(k)  [9], 

d  <9  d  E 

(2.25)  T(k)  =  - ci  ~k2^/kE(k)E(k)  +  c2  ^ k3  \/ kE(k)  — 

with  constants  compatible  with  the  existence  of  Kolmogorov  and  equipartition  spectra  [9].  The  evaluation 
of  the  moment  of  order  -1  has  been  proposed  [5]  for  the  EDQNM  energy  transfer  model. 

To  conclude,  write  the  e  transport  equation  in  the  form 

(2.26)  e  =  Ctllp-(CT-l)j 

The  values  of  C€\  and  Ct  are  subjects  for  future  investigation. 


2.3.  The  moment  equation  of  order  +2.  It  is  natural  to  attempt  to  derive  an  equation  for  viscous 
dissipation  by  multiplying  Eq.  (2.1)  by  vk 2  and  integrating.  As  noted  in  [6],  this  approach  has  the  very 
attractive  feature  that  the  divergences  [10]  of  order  Re1/2  cancel,  for 

r  OO 

(2.27)  /  k2T(k)  =  en2d 

Jo 


and 

(2.28)  —  f  dK  k2D(k)  =  —I'e2^3 k1^3 

Jo 

and  the  definition  of  Kd  shows  that  the  right  sides  of  Eqs.  (2.27)  and  (2.28)  can  sum  to  zero.  However, 


fOO 

(2.29)  /  dn  uk2U(k)  ~  vel/3S2K2J3 

Jo 

is  of  order  Re-1/2  and  the  remaining  contribution  from  the  transfer  term,  —vK^e  is  of  order  Re~l.  Thus, 
unless  other  contributions  can  be  found,  the  right  side  of  the  moment  equation  of  order  two  vanishes  in  the 
limit  of  infinite  Re.  Thus,  although  the  moment  equation  of  order  -1  does  not  provide  satisfactory  model 
constants,  it  generates  a  more  satisfactory  model  form  than  does  the  moment  equation  of  order  +2. 
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3.  Multiple-scale  Models.  The  simplest  generalization  of  the  previous  single-scale  model  results  from 
from  taking  the  piecewise-Kolmogorov  spectrum 

(3.1)  Ei{n)  =  e2/3K“5/3  for  K._1  <  k  <  Ki  with  i  >  1 


and 


(3.2)  T(k)  =  -  ^  ej[5(K  -  Ki-i)  -  5{k  -  «*)] 

l<t<n 

Eq.  (3.2)  defines  a  shell  model  of  turbulence,  in  which  energy  is  transferred  from  each  discrete  region  of  scales 
to  the  adjacent  region  of  smaller  scales.  This  picture  greatly  oversimplifies  the  actual  energy  transfer  process, 
which  is  mediated  by  triad  interactions  and  which  permits  both  forward  and  backward  energy  transfer  and 
transfer  between  non-adjacent  regions.  A  more  accurate  description  is  given  by  two-point  closures  like  DIA 
and  EDQNM  (Laporta,  1995). 

Set  corresponding  to  Eq.  (2.17), 

(3.3)  Ili (/c)  =  ^-Cd£1/3k~7/3S2  for  i  <  k  < 

O 

Then  the  moments  of  order  zero,  obtained  by  integrating  over  the  regions  <  k  <  give  partial  energy 
equations 


(3.4) 
where 

(3.5) 


h>i  —  q — i  q  4"  Pi  Pi 


P,  =  cdi)'%-T  -  «rv*) 

Di  =  -  <£/’) 


and  by  definition,  eo  =  0.  Unlike  the  single-scale  energy  equation  Eq.  (2.4),  Eq.  (3.4)  depends  on  the  special 
hypotheses  made  in  Eqs.  (3.1)- (3. 2). 

Eq.  (3.1)  implies 


(3.6) 


h  =  lcKen^{3-*;2/3) 


from  which  differentiation  in  time  gives 

(3.7)  k  =  CKe~1/3(K^{3  -  K~2/3)ii  -  Cxe^K-^ki-i  +  CK^/3K~5/3k,i 


As  in  the  previous  section,  another  relation  is  required  to  eliminate  the  new  quantities  and  we  can  proceed 
either  by  following  Schiesel’s  approach  of  introducing  a  new  differential  equation  for  ki  or  by  forming  moments 
of  the  the  dynamic  equations  in  each  spectral  region. 

3.1.  Schiestel’s  method.  There  are  several  natural  choices  to  close  ki  in  terms  of  local  turbulence 
quantities,  each  of  which  will  lead  to  a  different  model.  One  natural  generalization  of  Eq.  (2.12)  to  multiple 
regions  is  to  set 

(3.8)  ki  =  7 


so  that  the  ratio  ki/t^i  is  determined  by  the  local  Kolmogorov  frequency. 
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Schiestel  ([3],  Eq.  (24))  suggests  instead  that  ki  be  determined  as  the  ratio  of  the  flux  into  the  region 
Ki  <  k  <  Ki+ 1  divided  by  the  local  energy  density  at  Ki.  In  the  present  notation,  this  leads  to 


(3.9) 


5/3 

=  7 - 


'Ei+1(Ki)  ' ckel 


2/3 


■i+ 1 


where  the  choice  of  Ei+ i,  the  energy  density  for  scales  greater  than  Ki  seems  to  be  the  most  natural. 
Substituting  Eq.  (3.9)  in  Eq.  (3.7)  leads  to 


(3.10) 


•  2  ki  .  2/3 

ki  =  3 175- 

2  ei+l 


which  can  be  rearranged  as  the  partial  e  equation 


(3.11) 


3  /r>  ^  ,  3  3  e/3  e? 

€<  =  2F(Pi  "  +  2^  k~  “  0 

it  n>l  id  rvi  id  £  '  /Vj 

4+1 
,2/3 

i — 1 
2/3  J 
s+i 


3  e; 


=  2EPl+27 


3  €i—\€i  3  6, 


-2El1  +  7- 


This  result  differs  slightly  from  Eq.  (25)  of  [3]  because  of  the  final  term  which  modifies  the  destruction  term. 
Pi  can  be  neglected  except  in  the  regions  with  small  i  and  Di  can  be  neglected  except  in  the  regions  with 
large  i\  in  the  intermediate  regions,  the  energy  content  is  determined  by  nonlinear  energy  transfer  alone.  It 
is  therefore  a  reasonable  approximation  to  set  Pj  —  0  for  i  >  2  and  Di  —  0  for  i  <  n.  Then  the  production  of 
ti  comes  primarily  from  the  term  proportional  to  6*_i,  the  flux  into  the  region  k  >  Ki  from  the  region  k  <  Ki. 
With  these  approximations,  Eq.  (3.11)  reduces  to  the  transfer  equations  of  the  HLS  model  of  Eq.(l.l)  except 
for  the  modification  noted  above  of  the  destruction  term. 


3.2.  The  method  of  moments.  As  in  the  derivation  of  the  single-scale  model,  it  can  be  objected  that 
the  approach  of  [3]  requires  some  arbitrary  choices,  like  Eq.  (3.9).  A  natural  generalization  of  the  method 
of  moments  is  to  form  the  partial  moment  equations 


,3.12)  ±  f  r  *  £Wi  -  r  *  W  +  r  &  HiW  .  r  dK 

1  ^  ^  J  Ki-i  ^  J  1  ^  J  Ki- 1 

With  the  piecewise  Kolmogorov  forms  Eqs.  (3.1),  (3.2),  and  (3.3),  this  gives 


(3.13) 


i\o4r  wy-*-*)]- 

e»~l  ~  f»  ,  1/3,  -7/3  -7/3xc2 


+  f  Cfcer  KT  -  «r/3)52  -  3 uCk(kY3  -  k]L\) 

Ki—  l  ( 


AW 

K 


where  again,  e0  —  0.  I11  integrating  the  delta  functions  in  T;(k),  the  limits  are  1  <  k  <  Ki. 

The  elimination  of  the  quantities  Ki  is  straightforward  but  leads  to  lengthy  expressions.  To  illustrate  the 
procedure,  consider  a  two-scale  model  with  partial  energies  &i,/c2  and  partial  transfers  61,62.  Assume  that 
€2  balances  the  viscous  dissipation  so  that  this  model  generalizes  the  usual  two-equation  model  by  allowing 
transfer  and  dissipation  to  be  unequal:  compare  [6]. 

The  partial  energy  equations  are 


(3.14) 


h  =  Pi  -  d  -  Di 
k'2  =  P‘2  “f~  61  —  D'2 
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where  in  the  second  equation,  integration  over  k  >  K\  causes  the  contributions  from  e2  to  cancel.  Introduce 
the  approximations  discussed  above, 


P*Pi 
P2  »0 
D  D‘2 

(3.15)  A  «  0 
and  to  insure  the  overall  energy  balance 

(3.16)  k )1+k2  =  k  =  P-D 
set 


(3.17)  D  =  e2 
Then  Eq.  (3.14)  becomes 

hi  —  P  —  ti 

(3.18)  k2  =  e i  -  e2 


Effectively,  this  is  a  model  in  which  energy  transfer  through  the  large  scales  and  the  viscous  dissipation  can 
be  unequal.  However,  the  derivation  does  not  require  the  formation  of  the  moment  of  order  +2. 

Since 


(3.19) 


h  =  \cKeT^-K-^) 

k2  =  23Ki2/3 


the  scale  ratio  ki/kq  can  be  eliminated  through 


(3.20) 
or 

(3.21) 


Hsmr-1 


Note  that  Eq.  (3.20)  implies  that  k\  >  and  that  Eq.  (3.19  implies  that  Ki  =  k0  is  equivalent  to  ki  =  0. 
The  relations 

(3.22)  ki  =  -  Ck4/3ko 5/3*o  +  CK^3  *n>lz  ki 

o  €i 

(3.23)  k2  =  \  —e2  - 

o  e2 

can  be  used  to  eliminate  k0  and  k± . 

It  is  easiest  to  begin  the  derivation  with  the  second  wavenumber  partition.  A  calculation  similar  to  that 
leading  to  the  single-scale  result  Eq.  (2.26)  results  in 


(3.24) 


—  C2  +  k2  =  (ei  —  €2)Cti 
5  e2 
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where  the  present  theory  actually  predicts  CVi  =  1.  Even  if  this  value  is  left  general,  in  the  notation  of  Eq. 
(1.1),  this  calculation  predicts  Ct  1  =  C*2.  Eq.  (3.24)  can  be  rearranged  as 

(3.25)  e2  =  ~(1  -  Cti)(€i  -  e2) 

where  a  small  production  term  proportional  to  P2  has  been  neglected.  Note  that  the  production  of  e2  comes 
primarily  from  e\. 

Straightforward  calculation  leads  to  the  ei  equation 

(2ki  4/  k0\  k2 

I  5  ei  15  V  Ki  /  ej/3e^/3  / €l 

(3.26)  ~P  -  (1  -  Cro)ti  -  g  (l  -  )  (1  “  -  €2) 

Note  that  the  last  contribution  arises  from  the  k\  term  of  Eq.  (3.22),  which  is  eliminated  in  terms  of  e2 
through  Eq.  (3.23). 

To  summarize,  the  two-scale  model  contains  the  partial  energy  equations, 

*1  =  P  -  ci 

(3.27)  k2  =  ei  -  e2 


and  the  partial  dissipation  equations, 


f2fci  4/  kq\  k2  1. 

l5ei  15  V  kJ  z\/\y3  J£l_ 

(3.28)  ^  -  Cto^€i  ~  ^(1  _  (^)  (!  “  2^ri)(ei  -  €2) 


(3.29)  e2  =  ^g(l-C'T1)(e1-e2) 
with  the  definition 

(3.30)  l  =  {|(^/.  +  ir 

This  two-scale,  four-equation  model  separates  large-  and  small-scale  transfer  without  requiring  Reynolds 
number  dependent  coefficients.  Deficiencies  of  the  derivation  include  the  failure  to  predict  Ce2  and  the 
equality  Ct  1  =  Ct2  in  the  small-scale  transfer  equation  Eq.  (3.29). 

4.  Relaxation  to  the  Single-scale  Model.  Similarity  solutions  have  proven  indispensible  in  cali¬ 
brating  turbulence  models.  Examples  include  decaying  turbulence  and  spatially  self-similar  turbulent  shear 
flows  like  jets,  mixing  layers,  and  wakes.  All  of  these  flows  were  applied  to  calibrate  and  validate  the  first 
mixing-length  models  of  turbulence. 

We  must  expect  that  when  applied  to  a  turbulent  flow  which  relaxes  to  a  self-similar  evolution,  any 
multiple-scale  model  must  relax  to  a  single-scale  model  with  a  unique  similarity  solution.  In  particular,  the 
additional  freedom  allowed  in  multiple-scale  models  must  not  permit  spurious  results  like  multiple  power-law 
exponents  in  decaying  turbulence.  Both  the  HLS  model  and  the  proposed  model  will  next  be  analyzed  in 
decaying  turbulence,  and  the  conditions  which  prevent  the  existence  of  multiple  exponents  are  found.  In 
both  cases,  the  condition  is  a  simple  inequality  among  the  model  constants. 
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4.1.  The  HLS  model.  Consider  decaying  turbulence  described  by  the  HLS  model  Eq.  (1.1)  with 
P  =  0.  Look  for  a  solution 


(4.1) 


ft/p  — -  CLp  t 

kt  =  atta 
ep  -  epta~l 
et  =  ett a“1 


Substituting  in  Eq.  (1.1)  leads  to  a  system  of  homogeneous  equations  in  the  constants  ap)at,ep,et  which 
has  a  nontrivial  solution  provided 


(4.2) 


found  to  be 


a 

0 

i 

0 

0  = 

0 

a 

-i 

1 

a  -  1 

0 

C'P2 

0 

0 

a  -  1 

~Cn 

Ca 

-1/(^2 -1) 

and  a 

=  1/ {Cp2 

=  {aCt2  -  ol  +  1  ){aCP2(i  +  1) 


(4.3)  ap:  at- 

dp  I  (If  — 

The  first  solution  in  Eq.  (4.3)  obviously  corresponds  to  the  reduction  of  the  two-scale  model  to  a 
single-scale  model  since  kp  =  0.  To  avoid  the  existence  of  a  second  power  law  in  decaying  turbulence,  the 
second  solution  in  Eq.  (4.3)  must  be  non-realizable  or  unstable;  its  non-realizability  is  assured  if  one  of  the 
amplitudes  ap  or  at  must  be  negative.  This  occurs  if 


0:1  if  a  =  —  l/(Cj2  —  1) 

Cp2  —  Ct-2  ■  Ct2  —  Cn  if  a  =  —l/(Cp2  —  1) 


(4.4) 


Cp2  <  Cn  and  Cti  Cp2 


The  decay  equations  were  integrated  for  models  satisfying  and  violating  the  condition  Eq.  (4.4).  The 
results  are  shown  in  Fig.  (4.1).  First,  the  model  constants  were  arbitrarily  chosen  as  Cp2  =  1.5 ,Cn  = 
1-2, Ct2  =  2.0  and  the  initial  conditions  were  kp(0)  =  «t(0)  =  0.1,ep(0)  =  C((0)  =  1.0.  In  this  case,  Eq. 
(4.4)  is  satisfied.  The  resulting  decay  is  shown  in  the  left  graph  in  Fig.  (4.1).  It  shows  that  the  energy  kp 
approaches  zero  after  an  initial  transient,  indicating  that  at  long  times,  the  multiple-scale  model  reduces  to 
a  single-scale  model.  The  dotted  line  shows  the  power  law  decay  k,  ~  expected  in  this  case. 

To  demonstrate  that  Eq.  (4.4)  is  needed  because  the  second  solution  in  Eq.  (4.3)  can  be  stable,  the  decay 
equations  were  integrated  for  a  case  which  does  not  satisfy  the  constraint,  Cp2  =  3.0 ,Cn  =  1.2,  Ct2  =  2.0. 
The  results  are  shown  in  the  center  and  right-hand  graphs  in  Fig.  (4.1).  In  the  center  graph,  the  initial 
conditions  were  «p(0)  =  0,Kt(0)  =  0.1,ep(0)  =  0,et(0)  =  1.0  while  in  the  right-hand  graph,  the  initial 
conditions  were  np( 0)  =  ki(0)  =  0.1,ep(0)  =  et(0)  =  1.0.  Clearly,  the  power  law  decay  is  different  in  each 
case;  the  dotted  line  again  corresponds  to  the  power  law  k  ~  f — 1 .  Thus,  if  the  conditions  expressed  by  Eq. 
(4.4)  are  not  satisfied,  two  distinct  time-scaling  laws  can  exist  for  decaying  turbulence. 

4.2.  The  moment  model.  The  analysis  for  the  moment  model  Eqs.  (3.27)-(3.30)  is  similar.  It  is 
simpler  to  integrate  Eq.  (3.13)  including  the  empirical  factors  in  the  dissipation  rate  terms.  Substituting 
the  power-law  decay  Eq.  (4.1)  and  the  additional  relations 


(4.5) 


Ko  =  b0t  1 

Kl  =  blra>'2-1 
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turbulence  quantities 


Fig.  4.1.  decaying  turbulence ,  HLS  model:  (left)  Ct\  <  Cp 2  <  Ct2  one  decay  mode ;  (center)  Cn  <  Cp 2  <  Cf2  two  decay 
modes;  (right)  Cti  <  Ct2  <  Cp 2  two  decay  modes .  The  center  and  right  figures  show  the  possibility  of  two  different  power  law 
decay  rates  when  the  condition  in  Eq.  (4-4  n°t  satisfied. 


in  Eq.  (3.13), 

(4.6)  \cK^K,-'-p  -  i.rs/,)(|a  +  1)  =  -~Cri 

(4.7)  5ck4/V/s(5«+1)  =  -^ct! 


The  consistency  of  these  equations  evidently  requires 

(4.8)  k5/3  -  1  =  ^P-k(k2^3  -  1) 

Ct2 

where  k  —  ki/kq.  The  only  solution  is  k  =  1  corresponding  to  provided  Cti  <  Ct2  but  multiple 

solutions  are  possible  once  Cti  is  large  enough  relative  to  Ct2-  For  example,  a  solution  k  «  1.2  exists  if 
Ct2  =  2Cti-  Note  that  if  k  =  1,  then  k\  =  0,  again  indicating  reduction  of  the  multiple-scale  model  to  a 
single-scale  model. 

In  the  case  that  only  one  solution  exists,  the  power-law  decay  rate  is  again 


(4.9) 
where 

(4.10) 


1 


a  =  — 


<7,2-1 


C«2  =  jj(l-CTl) 


indicating  the  reduction  for  decaying  turbulence  to  the  correct  limit. 

Examples  of  decaying  turbulence  computed  with  the  moment  model  appear  in  Fig.  (4.2).  The  initial 
conditions  are 


left  graph  Fig.  (4.2):  Aq(0)  =  0.1  Jfe2(0)  =  0.1  ci(0)  =  1.0  c2(0)  =  1.0 

right  graph  Fig.  (4.2):  *i(0)  =  0.15  &2(0)  =  0.05  ci(0)  =  3.0  c2(0)  =  1.0 


These  initial  conditions  are  chosen  to  correspond  to  the  same  initial  conditions  in  a  single-scale  model, 
namely  k( 0)  =  fci(0)  +  ^(0)  —  0.2,  e(0)  =  e2(0)  =  1.0.  At  large  times,  both  calculations  follow  the  same 
power  law  decay,  but  there  are  clearly  differences  in  the  transient  evolution  before  self-similarity  is  obtained. 
The  decay  of  the  total  kinetic  energy  for  the  two  sets  of  initial  conditions  is  compared  in  Fig.  (4.3). 


ll 


Fig.  4.2.  decaying  turbulence,  moment  model :  (left)  initial  conditions  given  by  first  of  Eq.  (4>  11);  (right)  initial  conditions 
given  by  second  of  Eq.  (4-11).  The  graphs  show  the  effects  of  different  initial  conditions  on  the  decay  of  isotropic  turbulence 
predicted  by  a  two-scale  model.  The  different  sets  of  initial  conditions  correspond  to  the  same  initial  conditions  for  a  single-scale 
model 


Although  both  cases  follow  the  same  power  law  at  long  times,  transient  effects  cause  the  energy  to  decay 
more  quickly  for  the  second  set  of  initial  conditions  than  for  the  first  set.  The  capability  to  model  this  type 
of  transient  behavior  illustrates  the  justification  of  multiple-scale  modeling:  a  single-scale  model  could  not 
distinguish  between  these  two  cases. 


FlG.  4.3.  decaying  turbulence,  moment  model,  energy  decay  for  the  two  sets  of  initial  conditions  in  Eq.  (4-ll)> 

5.  Conclusions. 

1.  The  arbitrary  elements  of  the  derivation  of  multiple  scale  models  following  [3]  can  be  avoided  by  the 
method  of  moments,  which  reduces  the  continuous  spectral  evolution  equations  to  a  finite  dimen¬ 
sional  system. 

2.  The  moment  equation  of  order  -1  leads  to  a  well  defined  result,  but  the  moment  equation  of  order 
+2  is  problematical. 

3.  However,  the  moment  method  leads  to  difficulties  with  the  destruction  of  dissipation  term  both  in 
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the  single-scale  model  and  in  the  multiple-scale  models.  The  difficulty  originates  in  the  moment  of 
the  transfer  integral. 

4.  A  two  scale  model  in  which  energy  transfer  through  the  large  scales  can  be  distinguished  from  viscous 
dissipation  can  be  derived  without  forming  the  problematical  moment  of  order  +2. 

5.  Multiple-scale  models  must  reduce  to  single-scale  models  when  the  turbulence  evolution  is  self¬ 
similar.  Conditions  which  insure  this  reduction  are  derived  in  the  special  case  of  decaying  turbulence. 
The  same  analysis  should  be  completed  for  other  self-similar  flows. 
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